On the size scaling of the nearest level spacing at criticality 
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It is conjectured that the size scaling of the nearest level spacing in the critical spectral region, 
S{N) oc N~'^, remains qualitatively the same within phases of extended and critical states. The 
exponent A is therefore identical to that for the bare level spacing (at zero disorder). Our calcula- 
tion of the scaling for the one- dimensional model with diagonal disorder and long-range power-like 
interaction confirms the conjecture. 
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Localization-delocalization transition (LDT) in disor- 
dered systems, predicted by Anderson for 3D in 1958 0, 
still remains a thoroughly investigated problem (see 
Refs. for an overview). A remarkable progress has 

been achieved in understanding of critical systems, in 
particular, the multifractal nature of critical wave func- 
tions has been discovered d H 0, 13). Statistics of 
energy levels within disordered models have also been in- 
tensively investigated during the recent years. It was 
demonstrated that statistical properties of level spac- 
ing are closely related to general localization proper- 
ties 

The spacing between two energy levels is determined 
by the overlap of the corresponding wave functions. 
Thus, weak overlap of localized states results in Poisson 
statistics of energy levels in the thermodynamic limit, 
which implies that the level spacing can be arbitrary 
small and does not depend on the system size. Contrary 
to that, strong overlap of the extended eigenstates re- 
sults in strong level repulsion. Therefore, Wigner-Dvson 
statistics describe well energy levels in good metal j3.ll2j. 
The level spacing for extended states depends on the sys- 
tem size. In this contribution we consider the level spac- 
ing in a critical spectral regions, i.e. in a region where 
the localization-delocalization transition takes place or at 
the mobility edge. 

Within the phase of the extended states, the exponent 
A of the size scaling of the level spacing, S{N) oc N~^, 
is expected to be independent of disorder, otherwise the 
statistics of the energy levels spacing could change with 
disorder, which would be in direct contradiction to the ro- 
bustness of the Wigner- Dyson statistics within the phase 
of the extended states j3, 111, 01 ■ criticality, e.g. 
close to the transition point or the mobility edge, the 
level statistics, although manifesting strong level repul- 
sion, is different from the Wigner-Dyson one More- 
over, reasonings based on the consideration of wave func- 
tion overlap should be used with caution, because criti- 
cal eigenstates are very sparse and it is not obvious how 
much they overlap. Fyodorov and Mirlin addressed this 
problem within the framework of the power law random- 
banded model (see Refs. and references therein for 
the detailed description of the model). They showed an- 



alytically that two close critical eigenstates are strongly 
correlated, in the sense that the overlap of the two states 
is of the order of the self-overlap of a state 0|. It 
was conjectured in Ref. [l^j that the strong correlation 
of critical states is a general property. This correla- 
tion results in strong level repulsion up to the critical 
point [mil [13, El 111 Usui- Because the level spac- 
ing is determined by the overlap of wave functions, we 
conjecture that the strong overlap of critical eigenfunc- 
tions preserves the size scaling of the level spacing within 
the phase of the critical states also. Thus, the scaling ex- 
ponent A is expected to be the same throughout both 
phases of extended and critical states. The latter phase 
extends up to the magnitude of disorder at which the lo- 
calization length becomes equal to the linear system size. 
The exponent is therefore identical to that for the bare 
level spacing (at zero disorder) in the considered energy 
region. 

To prove the conjecture on the independence of the 
scaling exponent of disorder we consider the level spacing 
within the framework of the tight-binding Hamiltonian 
on a one-dimensional (ID) regular lattice of N sites with 
power-like inter-site interaction: 
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where \n) is the ket vector of a state with on-site energy 
En- The energies are stochastic variables, uncorrelated 
for different sites and distributed uniformly around zero 
within the interval of width W. The hopping integrals 
are J,nn — J/\m — n\'^, J„„ — with 1 < ^ < 3/2. We 
set J > 0, then the LDT occurs at the upper band edge 
provided l</i<3/2j2l[23. The latter fact makes the 
model very advantageous from the viewpoint of numeri- 
cal calculations because efficient algorithms, such as the 
Lanczos, can be used to find a few extreme eigenstates of 
interest. In all calculations we consider three uppermost 
states under open boundary conditions for system sizes 
in the range N = 4096 ^ 65536. All results are averaged 
over more than 510"^ x (65536/iV) disorder realizations. 
The lattice constant is set to unity. 

Recently, the model Q was studied numerically for 
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two values of interaction exponent /i: /i = 3/2 and 
/i = 4/3 |2J|. In the former marginal case no signs of 
transition was found, while at /i = 4/3 the model was 
found to reveal the LDT at = 10.9 ± 0.2. In this con- 
tribution we focus on the analysis of the level spacing. In 
particular, we study the size scaling of the level spacing 
at the top of the band in the vicinity of the critical value 
of disorder Wc for /z = 4/3. 

We are particularly interested in the scaling of the level 
spacing within the phase of the critical states. Because 
we deal with finite systems, the phase of critical states 
is of finite width, so we determine the width first. The 
critical phase extends from the critical point W = Wc up 
to the disorder magnitude W = Wl at which the local- 
ization length becomes equal to the linear system size: 
= L (within the considered ID model L = N). 
To determine the disorder magnitude Wl we analyze the 
size scaling of the typical value of the inverse participa- 
tion ratio (IPR). The IPR for the state a is defined in a 
standard way: 

N 
n=l 

where ipan is the a-th normalized eigenstate of the Hamil- 
tonian The typical value of the IPR is defined as 
I2 — cxp(ln/2"^), where the angle brackets stand for 
averaging over disorder realizations and a. We use the 
typical value of the IPR rather than the mean value of 
the IPR because, unlike the latter, the former is a self- 
averaged quantity 0, [2^ . 

Within the phase of extended states the IPR size scales 
as I2 oc , while for the localized states the scaling 
reads: I2 oc Wegner found that at criticality I2 size 
scales anomaluosly jj]: /20c L~^^ where the fractal di- 
mension D2 is smaller than the spacial dimensionality, 
which results from the multifractal structure of the eigen- 
functions 0, 0, IE • Above the transition point the 
scaling of the IPR is governed by the localization length ^ 
rather than by the system size: I2 oc . The local- 

ization length depends on the magnitude of disorder in 
the vicinity of the critical point ^{W) = io\W/Wc - l]^" 
{v > 0), therefore, two scaling regimes are expected 
above the critical point. Within the phase of the crit- 
ical states, i.e. for W > Wc and system sizes L < £,{W), 
the IPR scaling is still determined by the size of the sys- 
tem: I2 oc L~^^. Contrary to that, within the phase of 
the localized states, i.e. for L > £,{W), the IPR scaling 
is determined by the localization length £,{W) and does 
not depend on the system size. ^ 

Fig. n shows the size scaling of the quantity I2 — 
hL^'- The value of the correlation dimension D2 was 
obtained by studying the size scaling of the typical value 
of the IPR at criticality: D2 — 0.66 for /i = 4/3 (a 
detailed analysis of the correlation dimension is to be 
published elsewhere). Within the phase of the extended 
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FIG. 1: Size scaling of the parameter I2 — I2 in the 
vicinity of the critical point. 

states the quantity I2 decreases with the system size: 
I2 oc L^2-i < !)■ At the transition point, I2 

is expected to be constant. Above the transition, two 
scaling regimes are expected: I2 remains constant for 
L < £,{W), while for systems that are larger than the lo- 
calization length, I2 is expected to increase with the size: 
I2 oc ^ q^]^g above transparent scenario suggests 

that for a given set of system sizes the phase of the critical 
states extends from the critical point W = Wc up to the 
disorder magnitude W = W > Wc at which the size scal- 
ing of I2 starts to deviate from a constant (as a matter of 
fact, W = Wl for L = Lmax)- Fig. ^ demonstrates that, 
for the considered system sizes {L = N — 4096 ~ 65536), 
the critical phase extends up to the W « 11.5 J. 

Having obtained the disorder at which the uppermost 
eigenstates are critical we turn to the analyses of the size 
scaling of the level spacing for these states. Under peri- 
odic boundary conditions the level spacing at the upper 
edge of the bare band (for non-disordered system) de- 
pends on the system size N {N 3> 1) as follows: 

«.('V) = §^ + 0(^) . (3) 

where C{fi) = 2T{2 — fi) cos{tt {fi — l)/2). Under open 
boundary conditions, which we used in all calculations, 
no analytical expression for the spacing at the upper band 
edge can be obtained. Nevertheless, its size scaling is 
close to Eq. and the leading (non-zero) power of N 
in the expansion is also 1 — /.t. In particular, for fi — 4/3 
the spacing reads: 
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the constant C4/3 was obtained numerically: C4/3 « 3.58 
We calculated the size scaling of the mean level spacing 
at the upper band edge for fi = 4/3 and different values 
of disorder magnitude W. For each value of W, we fit the 
formula C TV"'*' to the calculated size dependence of the 
mean spacing, varying the system size within the range 
N = 4096 65536. The result of the fit is presented 
in Fig. 121 squares show the disorder dependence of the 
exponent A normalized to its value for non-disordered 
system, 1/3, while triangles show the disorder scaling of 
the factor C normalized to C4/3 (see Eq.^J. Vertical lines 
demonstrate the critical disorder VF=10.7-^11.1J [2^ . 
Error bars represent confidence intervals of the parame- 
ters C and A. The statistical error is of the order of the 
symbol size. 
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FIG. 2: Disorder scaling of the parameters A and C for fi — 
4/3 in the vicinity of the critical point. 

The figure demonstrates that within the phase of the 
extended and critical states, that is for VF < 11.5 J, 
the exponent that determines the size dependence of the 
mean level spacing in the critical energy region (i.e. at 
the upper band edge) is the same as that for the bare 
band: A = 1/3. The disorder reduces the level spacing, as 
expected, but does not affect the qualitative dependence 
of the spacing on the system size. Within the phase of 
the localized states, that is for W > 11.5 J, the local- 
ization length becomes smaller than the system size, and 
the level spacing is not determined by the size of the sys- 
tem. Therefore, the formula C N^^ becomes irrelevant 
within the phase of the localized states; the increase of 
the confidence intervals for the parameters C and A in 
the range W > 11.5 J manifests the irrelevance. 

To conclude, we have demonstrated that within both 
phases of the extended and critical states the size scaling 



of the level spacing remains qualitatively unchanged: the 
exponent of the size scaling does not depend on disorder. 
The exponent is therefore the same as that in the consid- 
ered energy region of the bare band (for a non-disordered 
system). The above result is probably general as it is a 
consequence of strong overlap of critical eigenstates that 
was reported on recently [ig. The independence of the 
scaling exponent of the disorder has several important 
consequences. It suggests, in particular, that the bare 
scaling of the energy levels can be used for qualitative 
argumentations that are applied to disordered systems 
(see, for example, Ref . and references therein). Prop- 
erties of critical states are usually studied by considering 
all states within an energy window in the vicinity of the 
critical point. The mean level spacing A oc A^^^ is of- 
ten used to estimate the number of critical states, which 
can lead to an over-estimation of this number, because 
the spacing in the critical region is much greater than A. 
In particular, for the standard three-dimensional Ander- 
son model the critical spectral region is the center of the 
band, where the spacing scales as 5* cx N^^^"^: the rele- 
vant spacing is therefore much greater than the average 
one. The over-estimation of the number of the critical 
states can lead to consideration of not only critical states 
but localized states also, which would result in the under- 
estimation of the correlation dimension, affect the energy 
level statistics, etc. 
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